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I .  INTRODUCTION 


The  purpose  of  this  report  is  to  document  the  significant 
results  of  the  second  year  of  our  research  program  under  Contract 
No.  DAAD05- 76-C-0757 .  The  object  of  this  portion  of  the  program 
was  to  apply  the  A.R.A.P.  Integral  Theory  of  Impact  to  the 
problem  of  modeling  the  behavior  of  long  rod  penetrators. 

4/ 

Over  the  past  several  years  the  Integral  Theory  of  Impact 
has  been  developed  by  A.R.A.P.  It  is  already  proving  to  be  a 

useful  tool  in  the  design  of  armor  and  penetrators.  The  theory 
contains  the  essential  physics  of  the  impact  process,  satisfies 
all  the  global  conservation  equations  and  is  contained  in  a 
computer  code  which  is  simple  and  inexpensive  to  operate.  The 
integral  theory  requires  far  less  empirical  information  than 
some  models  and  avoids  the  high  cost  and  complexity  of  multi¬ 
element  codes.  Its  simplicity  introduces  a  degree  of  economy 
that  makes  it  reasonable  to  conduct  parametric  studies  so  that 
predicted  trends  are  available,  rather  than  single  point  predic¬ 
tions.  This  aspect  of  the  simple  theory  greatly  facilitates 
the  interpretation  of  observations  and  the  selection  of  effective 
designs.  The  integral  theory  can,  therefore,  be  used  to  guide 
experimental  programs  and  to  select  those  designs  which  warrant 


further  study  using  the  large  codes. 

The  Integral  Theory  of  Impact  assumes  that  the  energy  of 
an  incoming  projectile  is  absorbed  by  xthe  armor  in  one  of  three 


forms : 


1.)  As  dissipated  energy  in  the  form  of  plastic  work  as 
the  target  flows  around  the  penetrator,  or  as  fracture  energy 
in  newly  created  fracture  surfaces .  This  nonrecoverable  portion 
of  the  total  energy  is  characterized  by  an  energy  per  unit  mass 
of  target  material,  E^.p  ,  and  it  appears  to  be  roughly  in¬ 
dependent  of  velocity  for  each  material.  The  product  of  target 
density  times  E*p  corresponds  to  the  "adiabatic  hardness"  of 
the  material,  or  its  hardness  measured  at  the  strain  rates  of 
impact . 


2. )  As  elastic  energy  absorbed  by  the  target  in  its  local 
elastic  deformation  near  the  penetrator  and  in  large  scale 
elastic  deformation  modes.  This  elastic,  or  recoverable,  por¬ 
tion  of  the  total  energy  is  parameterized  by  E*e  ,  the  elastic 
energy  per  unit  target  mass.  E,,rg  is  a  well  defined  function 
of  p/d  so  a  single  constant  defines  E^e  for  each  material 
over  the  complete  velocity  range. 

3. )  As  kinetic  energy  in  the  target  material  as  it  ac¬ 
celerates  and  begins  to  flow  around  the  incoming  penetrator. 

This  portion  of  the  energy  is  expressed  by  ,  where  V 

is  the  velocity  of  the  penetrator  face  relative  to  the  target 
and  Cj)  is  a  drag  coefficient  approximately  equal  to  1.0  for 
a  penetrator  with  a  spherical  front  end. 

Similarly,  when  analyzing  the  dynamics  of  the  deforming 
penetrator  itself,  there  will  be  a  quantity  which  measures  the 
dissipated  energy  per  unit  mass  absorbed  by  the  penetrator  as 
it  deforms  plastically  or  fractures.  It  is  the  analogue  of 
E*p  for  the  target  material.  We  will  call  this  quantity  E*^  . 

For  ductile  materials,  the  product  of  penetrator  density  times 
E*^  corresponds  to  the  "adiabatic  yield  strength"  of  the  material, 
which  is  the  uniaxial  yield  strength  of  the  material  measured  at 
the  strain  rates  of  impact.  E*^  is  assumed  to  be  a  constant 
for  each  penetrator  material. 

The  kinetic  energy  in  a  deforming  penetrator  is  modeled 
assuming  a  simple,  usually  linear,  flowfield  in  the  penetrator, 
and  simple  shapes  such  as  cubes  or  cylinders  to  approximate  the 
deformed  shape  of  the  projectile.  The  elastic  energy  in  pene- 
trators,  analogous  to  E*e  for  targets,  has  been  neglected  so 
far  because  it  is  relatively  small  compared  to  E^jj  at  the 
velocities  of  impact.  However,  in  principle  it  can  be  included 
also . 

Once  the  two  parameters  E^p  and  E^.e  are  known  for  a 
target  material,  and  E*^  is  known  for  the  penetrator,  the 
behavior  of  the  armor  and  penetrator  during  impact  can  be 


computed  from  global  energy  and  momentum  conservation  laws  by 
the  A.R.A.P.  Integral  Theory.  Since  it  is  only  the  sum  of  E*p 
and  E*e  which  governs  target  performance,  we  shall  often  refer 
to  the  sum  as  E,v  .  E*p  and  Ev.e  have  been  measured  in  im¬ 
pact  tests  for  a  variety  of  target  materials  from  lead  to  boron 
carbide,  from  salt  to  Rolled  Homogeneous  Armor,  over  a  velocity 
range  from  2b  ft/sec  to  6,400  ft/sec  and  have  been  shown  to 
provide  an  excellent  description  of  armor  behavior.^  Although 
E*  for  a  target  is  measured  in  impact  tests  with  nondeforming 
tungsten  carbide  balls,  the  same  value  of  E>v  for  the  target 
correctly  predict  its  performance  when  the  impactor  is  highly 
deforming,  such  as  lead  or  soft  aluminum  spheres^  or  with  a 
high  L/D,  such  as  a  long  rod  penetrator. 

Recently,  a  theory  has  been  developed  which  related  E*p 

and  E*e  to  more  fundamental  materials  properties,  such  as 

melting  temperature,  heat  capacity,  Young's  modulus  and  Brinell 
4 

hardness.  This  makes  possible  the  prediction  of  armor  performance 
from  static  tests  alone.  The  theory  has  been  verified  experi¬ 
mentally  over  the  same  wide  range  of  materials  for  which  impact 
experiments  have  been  conducted.  It  accurately  predicts  E,._ 
to  about  ±15  percent  for  all  these  materials.  This  formula 
has  enabled  A.R.A.P.  to  conduct  parametric  studies,  using  hand¬ 
book  properties  of  materials,  which  have  pointed  up  many  promising 
lightweight  armor  materials,  including  some  which  are  remarkably 
economical . 

In  this  paper,  we  shall  report  that  E^  for  penetrators 
can  also  be  computed  from  the  same  fundamental  materials 
properties,  and  can  therefore  be  predicted  from  purely  static 
tests.  Thus,  all  the  input  parameters  required  to  predict 
penetration  by  a  deforming  penetrator  into  a  target  can  be  ob¬ 
tained  from  handbook  values  of  materials  properties. 


II.  INTEGRAL  THEORY  FOR  ROD  PENETRATION 


The  Rod  Penetrator  Code  to  be  described  here,  which  we 
refer  to  as  "ROD,’’  uses  the  Integral  Theory  approach  outlined 
in  the  introduction.  As  a  penetrator  moves  through  a  target 
material  with  some  velocity  vface  at  the  penetrators  front  face, 
t he  Integral  Theory  for  target  performance  tells  us  that  the 
pressure  at  this  face  must  be 


Pressure 


-D  V2  +  F 
~T  ^  face  *t 


(1) 


where  is  the  target  density  and  PtL^t  the  adiabatic  hard¬ 
ness  of  the  target  material.  =  1  for  a  nondeforming 

spherical  front  face,  as  reported  previously.  This  formula  has 
been  verified  for  rigid  sphere  penetrators  and  deforming  sphere 
penetrators  over  a  wide  range  of  projectile  and  target  materials 
•uul  velocities.  Below  we  set  up  the  equations  governing  the 
internal  dynamics  of  a  long  rod  penetrator  which,  when  coupled 
with  Eq .  (1),  which  governs  the  target  dynamics,  completely 
specifies  the  problem. 

It  is  known  from  X-ray  photographs  that  the  stages  of  long 
rod  penetration  may  be  roughly  characterized  as  in  Fig.  1.  As 
t  tie  rod  impacts  the  target  the  pressures  generated  at  the  in¬ 
terface  begin  to  deform  the  front  end  of  the  penetrator,  is  in 
Fig.  lb.  Simultaneously  the  target  is  eroded  away  by  the  same 
pressure,  producing  a  crater,  as  in  Fig.  1c.  As  penetration  con¬ 
tinues,  material  at  the  leading  face  of  the  penetrator  is  eroded 
away  by  the  target,  forced  out  laterally  from  the  contact  region 
by  the  high  pressure  there  and  ejected  back  out  of  the  crater. 

As  material  is  eroded  from  the  rod  face,  new  material  is 
supplied  to  this  region  by  the  shaft  of  the  rod,  which  is 
traveling  at  a  higher  velocity  than  the  rod-target  interface. 

At  some  point,  Fig.  Id,  the  shaft  material  is  used  up  and  the 
head  is  decelerated  quickly  to  zero  velocity  by  the  target. 

We  shall  model  the  flowfield  of  a  rod  by  dividing  it  into 
two  regions;  the  head,  corresponding  to  the  front  region  of  a 
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Figure  2 


rod  which  is  undergoing  plastic  and  hydrodynamic  strain,  and 
a  rear  portion  of  the  rod,  the  shaft,  which  we  assume  to  be 
undeformed  rod  material.  In  Fig.  2,  the  assumptions  and 
parameters  of  the  model  are  summarized.  During  penetration 
the  head,  which  is  in  contact  with  the  target,  decelerates  and 
spreads  laterally.  We  assume  the  mass  flowfield  in  the  head 
is  linear,  and  the  head  is  of  cylindrical  shape.  The  motion 
of  the  material  in  the  head  is  characterized  by  a  center  of 
mass  velocity  ,  and  by  the  velocity  of  its  front  face  & 

and  side  face  b  relative  to  the  center  of  mass  of  the  head. 
Conservation  of  mass  across  the  boundary  between  the  shaft  and 
the  head  imposes  the  condition  that  the  rate  of  flow  of 
material  from  the  shaft  into  the  head  is 

=  rra2p  (V  -  Vf  +  Z)  (2) 

a  p  si 

where  a  is  the  radius  of  the  shaft,  p  the  penetrator  density, 

P 

and  Vg  the  velocity  of  the  shaft. 

As  penetration  proceeds,  the  head  widens  as  rod  material  is 
forced  to  flow  in  the  lateral  direction.  At  some  distance  from 
the  axis  of  the  rod,  say  cQa  ,  we  assume  the  flow  of  rod 
material  has  been  turned  or  sheared  off  by  the  target  and  no 
longer  can  apply  decelerating  forces  to  the  rod.  Thus,  when  rod 
material  in  the  head  moves  beyond  a  distance  e  a  ,  laterally, 
it  is  assumed  to  be  detached  from  the  rod.  The  dynamics  of  that 
material  as  it  is  further  slowed  by  the  target  will  not  affect 
the  deceleration  of  the  shaft  or  head.  This  assumption  is  justi¬ 
fied  for  some  cQ a  because  the  rod  material  at  this  point  in  the 
flow  has  been  adiabatically  heated  so  much  by  plastic  work  that 
its  shear  strength  is  very  low,  so  it  is  only  able  to  influence 
the  rod  through  compressive  or  hydrodynamic  forces.  However, 
the  axial  force  component  of  the  compressive  hydrodynamic  force 
on  the  rod  shaft  will  only  be  significant  within  one  or  two 
rod  radii  from  the  central  axis.  Thus,  we  shall  expect 
eoa  c  2a  .  Thus,  £Qa  really  characterizes  the  turning  radius 
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of  the  rod  material  in  the  target  or  the  shape  of  the  flowfield 
in  the  head.  We  further  assume  for  simplicity  that  the  shape  of 
the  flowfield  in  the  head  region  of  the  rod  does  not  change  too 
much  from  material  to  material.  Therefore,  we  may  take  eo  to 
be  the  same  constant  for  all  rod  penetrators,  no  matter  what 
material . 

When  b  (the  radius  of  the  rod  head)  reaches  t-.  a  ,  the 

o 

cut-off  radius,  we  assume  any  further  increase  in  the  radius  of 
the  head  simply  results  in  loss  of  rod  material  across  the 
boundary  at  eQa  ,  as  in  Fig.  2.  The  rate  of  mass  loss  from 
the  head  will  be 

Mb  =  47r£bppb  (3) 

where  b  is  the  lateral  velocity  of  material  in  the  head  of  the 
radius  b . .  Then  the  rate  of  change  of  mass  in  the  head  is 

3E  Mf  '  “a  -  Mb  « 


and  the  rate  of  change  of  mass  in  the  shaft  is 


s 


M 


a 


(5) 


The  pressure  applied  to  the  rod  front  face  by  the  target 
during  penetration  is 

C 

+  "2^  (Vf  +  &)  )  (6) 

as  described  in  Eq.  (1),  where  the  front  face  of  the  rod  moves 

at  the  velocity  +  £  .  This  pressure  acts  across  the  entire 

frontal  area  of  the  rod  head  which  is  in  contact  with  the 

2 

target.  The  area  is  rrb  giving 

C 

Total  force  =  ~  (MfVf  +  MgVs)  =  -nb2p  (E*t  +  ^  (Vf+;)2)  -  MbVf 

(7) 

from  total  conservation  of  momentum.  The  second  term  on  the  right 
accounts  for  momentum  loss  through  mass  loss  out  the  side  of  the 
head.  Substituting  from  (4),  we  find 

c 

Ma (Vf  -  Vs)  +  MfVf  +  MgVs  =  -rb2pt((E*t  +  ^(Vf+£)2)  (8) 
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This  may  be  separated  into  two  equations  for  the  conservation 
of  momentum  of  the  head  and  the  shaft.  The  shaft  will  only 
experience  decelerating  forces  if  it  has  a  nonzero  yield 
strength  a  .  In  this  case,  we  have 

(Force  on  Shaft)  =  -ira  o  =  M  V  .  (9) 

s  s 


When  o  =  0  ,  as  in  the  case  of  a  shaped  charge  jet,  which  is 
liquid,  the  shaft  velocity  remains  constant  throughout  penetra¬ 
tion.  Subtracting  this  from  (8)  above,  the  corresponding  equa¬ 
tion  for  the  head  is  obtained: 

O  C 

MfVf  =  -TTb2pt(E*t  +  2^(Vf+i)2)  +  OTia2  +  Ma(Vg  -  Vf)  (10) 


The  first  term  on  the  right  is  the  force  on  the  head  due  to  the 
target  pressure,  the  second  term  is  the  acceleration  of  the 
head  due  to  the  push  from  behind  applied  by  the  yield  strength 
of  the  shaft,  and  the  third  term  is  the  momentum  added  to  the 
head  from  the  material  passing  into  the  head  from  the  shaft. 


The  equations  above  account  for  momentum  conservation. 
Next,  we  require  energy  conservation.  The  total  kinetic  energy 
in  the  rod  is  given  by 


K 


Vl 

M  -# 
s  2 


+  l  b2) 


(ID 


where  the  third  term  on  the  right  accounts  for  the  internal 
kinetic  energy  in  the  flowfield  of  the  head,  a  =  3  for  a 
cylindrical  head.  The  total  work  per  unit  time  done  on  the 
rod  by  the  target  is 


U  =  Pt^b2(Vf+l)  (E*t  +  ^(Vf+£)2) 


(12) 


This  work  is  converted  into  either  heating  of  the  rod  of 
changing  the  kinetic  energy  of  the  rod.  The  heating  rate  is 
given  by  the  W  ,  the  rate  at  which  rod  material  is  converted 
into  the  hydrodynamic  state,  given  by 


(13) 
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where  is  the  "adiabatic  yield  strength"  of  the  rod 

material,  and  M  is  the  rate  at  which  rod  material  enters  the 

Si 

head.  The  conservation  of  total  energy  requires: 

W  +  K  +  U  =  0  (14) 

Lastly,  the  parameter  i  for  the  half-thickness  of  the  head  and 
b  for  the  radius  of  the  head  are  related  to  the  mass  of  the 
head  by 

27:ib^Pp  =  (15) 

These  equations  completely  specify  the  problem.  We  solve 
them  to  obtain  the  following  coupled  set  of  differential 
equations . 

M-  «  M  -  M,  (16) 

tab 

M  =  ira2p  (V  -  V-  +  i)  *  f,  (17) 

3  p  S  X  1 

Mb  =  4TTibppb  '  f2  (18) 

Q 

v£ '  -  fc(’b2pt(E-t  +  2s  <Vi)2)  +  Vvf  -  V  -f2  ■  °’a2) 


V  =  -a-rra  /M 
s  s 


C  .  i'l  O 

K  =  -(Vf+£)Trb2pt  (Vf+£) 2  4-  E*t)  -  MaP  E*d  -  -j-  (vf  +  b) 

/ w  X  (21) 


b  _  1  \  Mf  "  4 
b  =  K4  +  K5£ 


where 


-  (1  -  f2) 


2'  b 


b  b  f 
I  ’  1  M 


Mf\2 


f1T.a2o  <VS-Vf)  +  f24irlbp  b  (f) 
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and 


K 


5 


i 


where 


Q2 


a2fi  i  b 
Sib  1  i 

o  Mf  • 

K  +  Q2  - 1  r  bK4 

T-  +  7  bK5> 


and  penetration  rate 
P  =  Vf  +  i 


(25) 


(26) 


2 

M  V 
a  s 


^  (b2  +  V2)  +  MfVfVf| 


(27) 

(28) 


and  f2  are  integer  quantities  introduced  to  allow  the 
numerical  integration  to  proceed  smoothly  at  certain  discontinuous 
transition  points  in  the  model. 

f,  remains  1  unless  the  mass  M  of  the  shaft  becomes 
i  s 

zero.  This  will  happen  when  the  shaft  has  been  consumed  by 

erosion.  When  the  mass  M  becomes  zero,  f,  =  0  . 

s  1 

f2  remains  zero  until  b  ,  the  head  radius ,  reaches  eQa  . 
Then  f2  becomes  1  ,  and  prevents  the  radius  of  the  head  from 
increasing  beyond  eQa  .  Thus,  mass  loss  at  the  head  also 
begins  to  occur  when  f2  =  1  . 

This  set  of  equations  is  incorporated  into  the  computer  code 
ROD,  which  is  reproduced  in  Appendix  I.  The  input  parameters 
required  to  operate  the  code  consist  only  of  the  length  and 
radius  of  the  rod,  the  density  of  the  rod  material,  and  its 
adiabatic  hardness  pE*  ,  plus  corresponding  quantities  for  the 
target.  The  other  parameters  in  the  set  of  equations  above  are 
disposed  of  in  the  following  way:  We  have  learned  that  o  ,  the 
yield  strength  of  the  rod,  is  just  its  "adiabatic  uniaxial  flow 
stress"  or 


0  =  PpE*d 


(29) 


The  cut-off  radius  e  is  assumed  to  be  a  constant,  e  =  1.36  , 

o  o 

for  all  materials.  Furthermore,  the  penetration  depth  in  the 
rod  program  is  not  very  sensitive  to  the  initial  assumed  value 
of  i  ,  so  we  always  set  JL  =  a  ,  initially,  but  this  assumption 
is  not  critical.  The  only  input  parameters  needed  to  operate 
the  code  are  the  physical  dimensions ,  velocities ,  and  densities 
of  the  target  and  penetrators ,  plus  the  E,,.  values  of  the 
materials.  The  value  of  the  plastic  component  of  E,..  for  any 
material  can  be  derived  from  the  formula 


E„,.  =  0.55  C  T  £n 

pm 


°f(V£> 

•08pcptm 


as  derived  in  our  previous  interim  report,  where 
head  capacity,  T  the  melting  temperature,  p  t 


i  C  is  the 
P 

the  density 


head  capacity,  T  the  melting  temperature,  p  the  density 

and  a*  the  strain  rate  corrected  flow  stress  of  the  material. 
F 

Figure  3  displays  the  value  of  EVc  predicted  by  this  formula 
as  a  function  of  Brinell  hardness  for  a  number  of  materials 
of  interest.  When  elastic  effects  can  be  neglected  as  they  can 
be  for  most  armor  materials,  we  simply  take  E*t  of  the  target 
equal  to  E*  in  Eq.  (30). 

As  we  shall  show  in  this  report,  for  penetrators,  the 
corresponding  E*^  may  be  found  from 


E*d  “ 


where  x  =  0.42  for  the  code  ROD  and  E*  for  the  rod  is  com¬ 
puted  by  substituting  the  melting  temperature,  heat  capacity, 
density  and  flow  stress  values  of  the  rod  material  into  (30) . 


The  values  of  the  constants  in  (30)  have  been  modified  slightly 
from  those  given  in  Reference  1,  as  a  result  of  more  extensive 
impact  data. 


THEORETICAL  E*  vs  BRINNELL  HARDNESS 
FOR  ASSORTED  MATERIALS 


III.  DISCUSSION  OF  THE  SIMPLIFIED  CODE  "PEN" 


Recently  a  second  model  for  rod  penetration  has  been  developed 
which  is  conceptually  equivalent  to  the  code  ROD,  but  requires 
considerably  less  computer  running  time  and  is  more  useful  for 
analytical  investigations  because  of  the  simplicity  of  the  equa¬ 
tions.  Most  of  the  equations  (16)  -  (28)  are  used  in  calculating 
the  acceleration,  energy  and  momentum  in  the  flowfield  of  the  head. 
These  can  be  replaced  by  a  term  in  the  equation  of  momentum  con¬ 
servation  (19)  which  accounts  for  accelerations  in  the  head.  The 
resulting  model  and  code  is  named  "PEN." 

As  in  the  ROD  code,  a  two-element  rod  is  assumed,  consisting 
of  a  shaft  and  a  head.  The  cutoff  radius  eQa  is  retained,  and 
again  is  assumed  to  have  a  constant  value,  independent  of  rod  mat¬ 
erial.  The  initial  length  of  the  head  region  is  assumed  to  be 
a  ,  and  the  volume  of  the  head  is  assumed  to  be  constant  as  it 
flattens  and  widens  upon  impact.  Thus, 


i* 

“f 


ma 3  p. 


(32) 


The  decelerating  force  acting  on  the  head  is  the  sun  of  the 
decelerating  forces  provided  by  the  target  plus  the  accelerating 
force  on  the  head  supplied  by  the  yield  strength  of  the  shaft: 


MfVf  = 


-  mb 


!pt(E*t 


Vz  ) 

tV 


*a2pp(E*d  +  (Vg  -  Vf)2)  (33) 


where  the  first  term  on  the  right  is  the  force  exerted  on  the 
head  by  the  target  and  the  second  term  is  the  force  which  the 
shaft  exerts  on  the  head.  The  term  p  (V  -  Vf) 2  accounts  for 
the  momentum  gained  by  the  head  from  material  which  has  entered 
the  head  from  the  shaft.  The  front  face  velocity  of  the  rod  is 
assumed  to  be  equal  to  the  center  of  mass  velocity  of  the  head, 
.  Substituting  Eq .  (32)  into  Eq .  (33)  produces 
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Of  ■  iF-  (-»t  *7  CE*t  +  TVf>  +  cp  <<Vs  -  Vf>!  +  E*d» 

p  a  r 

(34) 

The  corresponding  equation  for  the  velocity  of  the  shaft  is 
Eq.  (20),  which  we  rewrite  as 

Vs  =  '  ^  na.2o  (35) 

where 

a  =  Pp  E*d  • 

The  penetration  p  of  the  rod  is  given  by 

P  -  Vf  ,  (36) 

and  the  erosion  of  the  rod  length  L  is  governed  by 

L  =  Vf  -  Vs  .  (37) 

The  radius  of  the  head  b  during  the  early  stages  of  pene¬ 
tration  when  it  is  widening,  is  computed  in  the  following  way: 

If  £  is  the  thickness  of  the  head, 


yielding 


frb2  £Pp  =  Mf  =  pp  7Ta; 


But  since  the  initial  thickness  of  the  head  is  a  ,  then  £  = 

a  -  (L  -  L)  ,  for  b  <  e  a  ,  so 
o  o 


)  (38) 
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defines  b  and  completes  the  system  of  equations.  This  system 

of  equations,  Eq .  (34)  -  Eq.  (38),  defines  the  numerical  code 

PEN.  As  with  the  ROD  code,  the  only  input  parameters  required 

are  the  initial  velocity  of  the  rod,  the  values  of  E,v  for 

the  rod  and  target,  and  the  dimensions  and  densities  of  the  materials 

This  model  can  be  related  to  hydrodynamic  models  of  rod  pene- 

trators  in  the  following  way.  During  steady-state  penetration 

of  a  rod,  typically  approaches  some  constant  value  ~  i  V  - 

.1  .  j  s 

j  Vg  .  During  this  stage  of  penetration,  may  be  neglected 

in  Eq.  (33).  Assuming  b  has  reached  its  full  value  of  eQa  , 

Eq.  (34)  may  be  approximated  as 

°D 

pP(Vs  ‘  Vf>*  =  ptciTVf  +  (pteoE*t  -  ppE*d>  (39) 

This  equation  may  be  compared  to  various  models  for  rod  pene¬ 
tration,  such  as  that  found  in  Ref.  6,7 


2  PP(VS  -  vf)>  =.  \  ptvj  +  <R  -  Y>  , 


where  R  is  the  target  strength  and  Y  the  strength  of  the  pene- 
trator.  Dividing  Eq .  (39)  by  2  and  comparing  coefficients,  we 
find 

1  ezCn  < - *  1 

2  o  D 

pteoE*t  //1N 

-  < - >  R  (41) 

PDE*d 

-£Y~  < - *  Y 

Now  the  values  of  R  and  Y  which  give  the  best  fit  to  experi¬ 
ment  for  a  number  of  materials  have  been  deduced  in  Refs.  (6-8). 

In  Table  1,  we  compile  the  experimental  values of  R  and  Y 
as  well  as  our  theoretical  prediction  of  them  based  on  Eq.  (41) 
and  the  theoretical  value  for  E*  ,  Eq .  (30).  The  good  agreement 
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indicated  by  Table  1  strongly  indicates  that  the  parameters  R 
and  Y  employed  in  the  hydrodynamic  models  are  related  to  the 
single  fundamental  material  property  E*  by  Eq .  (41 ) - 

Conceptually,  then,  the  rod  program  PEN  is  equivalent  to  a 
hydrodynamic  model  with  strength  in  the  target  and  rod,  such  as 
Eq .  (40),  plus  an  additional  term  proportional  to  the  decelera¬ 
tion  of  the  rod  front  face,  which  takes  into  account  inertial 
effects  at  the  front  of  the  rod  during  the  early  stages  of  impact 
before  equilibrium  of  pressures  has  been  established.  This  in- 
ertial  term  MfVf  acts  as  an  effective  stiffening  or  strengthen¬ 
ing  of  the  rod  during  this  early  phase  of  penetration,  and  accounts 
for  the  observed  fact  that  rods  of  lower  L/D  have  a  greater  pene¬ 
tration  vs.  rod  length,  P/L  ,  than  large  L/D  rods  do.  Were  it  not 
for  this  term,  and  an  additional  small  effect  due  to  the  initial 
spreading  of  front  face,  rods  of  all  L/D  would  have  the  same 
penetration  vs.  length  at  the  same  velocity,  as  the  hydrodynamic 
theory  of  Eq.  (40)  predicts. 

The  program  PEN  has  been  extremely  useful  as  a  means  of  gain¬ 
ing  an  intuitive  understanding  of  rod  penetration,  since  it  can 
predict  rod  performance  over  the  same  range  of  materials  and  vel¬ 
ocities  as  the  ROD  code,  yet  has  simpler  equations  which  can  be 
dealt  with  and  understood  algebraically. 

In  Section  VI,  the  predictions  of  PEN  are  compared  to  the 
code  ROD  and  to  experimental  data  for  a  wide  range  of  materials 

and  velocities.  We  have  found  that  the  value  of  e  and  v 

o  ^ 

which  gives  the  best  fit  to  experiment  for  PEN  is  e  «  1.7  , 

X  =  1.0  ,  and  =  0 . 5  .  These  values  are  used  in  all  computa¬ 
tions  employing  the  PEN  code. 

The  relationship  between  the  material  strengths  Y  and  R 
of  Table  1  is  reminiscent  of  the  relationship  between  the  uni¬ 
axial  tensile  strength  and  the  Brinell  hardness  of  a  material 
in  static  tests.  It  is  well  known  that  the  Brinell  hardness  B 
for  a  ductile  material,  which  is  just  the  pressure  that  a  ductile 
target  can  sustain  when  impressed  by  a  rigid  ball  indenter,  may 
be  related  to  the  uniaxial  flow  stress  of  a  rod  of  the  same  material 
by 
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■7  =  (0 . 3)  B  . 

Similarly,  the  effective  strength  R  of  the  target,  which 
is  analogous  to  B  ,  may  be  related  to  the  uniaxial  strength  Y 
of  the  rod,  which  is  analogous  to  o  ,  by 


.THEOR . 


E,./2 


^THEOR-  "  PpETTy 


,THEOR . 


=  .34  R 


THEOR. 


This  strongly  suggests  that  the  relationship  between  Y  and  R 
is  just  that  between  a  uniaxial  tensile  test  and  a  Brinell 
hardness  test  done  at  the  strain  rates  of  impact.  Thus  pE.  , 


which  determines  both 


is  a  true  measure  of  a 


materials'  strength  at  impact  strain  rates.  The  shear  heating 
process  in  the  deforming  material  at  these  strain  rates  is 
adiabatic  rather  than  isothermal,  since  heat  is  generated  locally 
in  regions  of  shear  much  faster  than  it  can  dissipate  by  thermal 
conduction.  When  the  local  heating,  with  the  attendant  local 
softening  of  the  material,  is  considered,  as  we  discussed  in 
Reference  4,  Eq .  (30)  is  derived  for  the  effective  material 
strength  at  impact  strain  rates  .  We  refer  to  pE.v  as  the 
"adiabatic  hardness"  of  a  material,  and  note  that  it 
determines  both  target  and  rod  strengths  at  these  strain  rates. 

The  quasi-hydrodynamic  model  of  Eq .  (39)  is  useful  in 
another  way  as  a  tool  for  analyzing  qualitatively  different 
regimes  of  penetration.  If  the  target  is  very  hard  (high  E,,.t  ) 
then  the  rod  will  not  penetrate  unless  its  velocity  is  suffi¬ 
ciently  high  to  overcome  the  target  strength  with  kinetic 
energy.  The  condition  for  the  lower  limit  of  velocity  required 
for  penetration  is  formed  by  setting  =  0  in  (39)  and 
solving  for  V  • 


(V0  - 
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Penetration  will  not  occur  unless  the  initial  rod  velocity 
exceeds  (Vs)q  >  according  to  this  model.  In  reality,  some 
penetration  does  occur  below  this  velocity  but,  as  in  Figs.  31, 

32  and  33,  there  usually  is  a  long  straight  section  of  the  P/L 
vs  velocity  curve  which,  when  extrapolated  to  zero  penetration, 
intersects  the  velocity  coordinate  at  a  value  given  approximately 
by  (Vs)  .  See  Ref.  6  for  a  discussion  of  this  relating  to 
the  data  in  Fig.  31.  The  discrepancy  at  velocities  below  (Vs) 
occurs  because  we  have  neglected  the  term  in  Eq .  (39). 

When  the  rod  is  very  strong  compared  to  the  target,  it 
may  not  erode  at  all,  and  then  it  behaves  like  a  nondeforming 
rod.  This  limit,  in  which  Vf  =  Vg  ,  will  occur  when 


If  the  quantity  under  the  radical  is  >0  ,  there  will 
exist  a  value  (Vs)  for  which  penetration  of  the  rod  can 
occur  without  erosion.  For  any  initial  velocity  below  (Vs)0 
the  rod  behaves  as  a  nondeforming  projectile,  and  the  A.R.A.P. 
integral  theory  for  nondeforming  projectiles  is  employed  instead 
of  the  deforming  rod  equation.  For  velocities  above  (Vs)q 
the  deforming  rod  model  applies.  Figure  35,  discussed  in 
Section  VI,  affords  a  striking  example  of  the  transition  from 
nondeforming  penetration  to  deforming  rod  penetration  as  the 
striking  velocity  is  increased. 


I 


IV .  BACKFACE  EFFECTS 

In  order  to  apply  the  A.R.A.P.  integral  theory  to  targets 
of  finite  thickness,  backface  effects  must  be  included.  The  E,.. 
concept  was  originally  developed  for  the  flow  of  target  material 
around  a  penetrator  in  a  semi-infinite  target.  The  shear  work 
done  on  the  target  material  in  the  flow  volume  defines  E,vt  , 
the  E*  for  the  target.  When  the  projectile  has  penetrated  al¬ 
most  all  the  way  through  the  target,  to  within  one  or  two  dia¬ 
meters  of  the  backface,  the  target  material  can  spall  or  simply 
bulge  on  the  backside,  rather  than  flowing  around  the  penetrator 
hydrodynamically .  Thus,  each  small  volume  of  target  material 
absorbs  less  energy  than  it  would  in  the  semi-infinite  case. 

Thus,  the  effective  E,v  for  the  target  decreases  near  the  back- 
face,  and  we  call  this  the  backface  effect. 

In  order  to  characterize  the  backface  effect  empirically, 
static  Brinell  hardness  tests  were  performed  on  1100-F  aluminum 
and  lexan  sheets  unsupported  at  the  back,  using  a  .250"  diameter 
WC  ball  at  very  shallow  and  very  deep  penetrations.  At  the  deep 
penetrations,  backface  effects  in  the  target  sample  affected  the 
hardness  measurement.  In  Fig.  4,  the  Brinell  hardness  of  a  1/4" 
thick  1100-F  aluminum  plate  is  measured  vs.  penetration  depth 
of  the  ball,  and  plotted  as  curve  (a).  The  hardness  is  roughly 
constant  with  penetration  until  the  ball  is  about  .150”  from  the 
backface,  at  which  point  the  hardness  begins  to  decrease  linearly 
with  further  penetration.  When  the  front  face  of  the  ball  reaches 
.475”  of  "penetration,"  so  it  has  actually  passed  through  the 
plate,  the  bulge  on  the  back  of  the  plate  fractures  and  the  mea¬ 
sured  hardness  drops  to  zero.  For  comparison,  the  hardness  vs. 
penetration  depth  for  the  same  WC  ball  in  a  semi-infinite  1100-F 
aluminum  plate  is  plotted  as  curve  (b) .  The  ratio  of  measured 
hardness  in  the  .250"  thick  plate  to  the  measured  hardness  in  the 
semi-infinite  plate  is  plotted  on  the  same  graph  as  a  solid  curve, 
(c) .  It  is  clear  that  as  the  ball  approaches  within  a  diameter 
or  so  cf  the  backface,  the  hardness  begins  to  decrease  mcnotoni- 
cally  with  penetration.  We  should  expect  that  pE*  ,  which 
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Finite  target  hardness 
Semi-infinite  target  hardness 


measures  the  target  strength  at  the  strain  rates  of  impact, 
should  decrease  in  roughly  the  same  way  as  the  measured  static 
hardness  does  near  the  backface. 


In  a  second  experiment,  the  Brinell  hardness  at  deer  pene¬ 
tration  in  lexar.  plates  was  measured  over  a  range  of  plate  thick¬ 
ness  and  ball  diameters.  The  results  are  shown  in  Fig.  5.  It 
was  found  that  the  results  could  be  fit  empirically  by 


Brinell  hardness  B  = 


f  Bo  ,  t  ■  p  >  (:  •  1): 


t  -  p  +  r 


Bo,i-p<(S-l)r 


where  £  =  A  ,  Bq  is  the  hardness  of  a  semi-infinite  lexan  plate, 
t  is  the  plate  thickness,  p  the  penetration  denth  of  the  ball, 
and  r  the  contact  radius  of  the  depression  made  by  the  ball  in 
the  target.  Obviously,  r  <.  a  ,  where  a  is  the  ball  radius,  and 


p  <  a 
p  >  a 


(A3) 


In  Fig.  6,  this  expression  is  compared  with  the  measured  value  of 
relative  hardness  at  various  penetrations  for  the  aluminum  plate 
discussed  in  Fig.  A.  The  agreement  is  qualitatively  good,  al¬ 
though  there  are  certainly  other  expressions  which  would  charac¬ 
terize  the  hardness  near  the  backface  as  well. 


The  form  of  expression,  (Eq .  (A2) ,  was  chosen  as  our  model 
for  the  backface  because  of  the  following  intuitive  model  for 
backface  effects.  We  assume  that  the  flowfield  of  target 
materials  around  the  penetrator  extends  for  some  distance  in  front 
of  the  penetrator.  We  expect  this  distance  to  be  proportional  to 
the  contact  radius  r  of  the  penetrator,  and  to  extend  a  distance 
£r  in  front  of  the  penetrator,  from  the  point  of  maximum  contact 
width.  Thus,  for  a  spherical  indenter  imbedded  less  then  one 
radius  deep  in  a  target,  the  flowfield  is  assumed  to  extend  a 
distance  gr  from  the  surface  of  the  target.  Once  the  ball 
is  fully  imbedded,  then  r  =  a  ,  and  the  flowfield  extends  a 
distance  ga  in  front  of  the  ball,  as  measured  from  the  center 
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Brinneil  hardness,  kg/mm 
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Figure  5 
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COMPARISON  OF  BACKFACE  MODEL  FOR  E#  WIT 
RELATIVE  HARDNESS  OF  ALUMINUM  MEAS- 


Figure 


Residual  velocity,  m/sec 


COMPARISON  OF  BACKFACE  MODEL 
WITH  EXPERIMENT 


Target  thickness,  inches 


Figure  7 


of  gravity  of  the  ball.  Thus,  the  backface  effect  will  begin  to 
show  up  when  p  -  a  =  t  -  6a  for  a  fully  imbedded  ball,  hence 
Eq.  (42).  For  a  blunt  nosed  object,  such  as  a  cube  or  cylinder 

impacting  end  on  into  a  target,  the  widest  point  of  the  penetra- 
tor  occurs  at  the  leading  face,  unlike  the  sphere.  Then  the 
backface  effect  begins  to  occur  when 


so  we  can  summarize : 


p  =  t  -  6a  , 


Backface  Nfodel 


For  a  sphere : 


E 


*t 


E*  ,  x  -  p  >  (6  -  1) r 

E*  (T  ,r )  ,  t  -  p  <  (6  -  1) r 


(44) 


For  a  cube  or  cylinder : 


,  t  -  p  >.  Br 
,  t  -  p  <  Sr 


(45) 


We  have  made  the  assumption  here  that  the  dynamic  stength  of  the 
target,  E*  ,  decreases  near  the  backface  in  the  same  way  as  the 
static  strength  does.  We  do  not  a  priori  expect  these  formulas 
with  the  same  constants  to  work  well  for  all  materials,  since 
brittle  materials  will  show  spall  and  other  backface  effects  when 
the  penetrator  is  many  diameters  from  the  backface.  It  is  pos¬ 
sible  that  by  making  6  inversely  proportional  to  the  failure 
strain  of  the  target,  the  formula  may  be  generalized.  Such  ap¬ 
proaches  will  be  considered  in  subsequent  work.  For  many  ductile 
materials,  however,  we  have  found  that  Eq.  (44)  and  Eq .  (45) 
accurately  describe  the  decrease  of  E*  near  the  backface. 

As  an  example,  refer  to  Figs.  7,  8  and  9  in  which  nonde¬ 
forming  tungsten  carbide  balls  and  highly  deforming  lead  projec¬ 
tiles  were  fired  into  rolled  homogeneous  armor,  5083  aluminum, 
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COMPARISON  OF  BACKFACE  MODEL 
WITH  EXPERIMENT 

WC  — -  ALUMINUM 
-  E  *t=  2  45  Btu  /  Ibm 


and  titanium  armor.  In  each  case  the  theoretical  residual  ve¬ 
locity  of  the  projectile  is  plotted  vs.  target  thickness  for  a 
given  initial  projectile  velocity,  using  the  A.R.A.P.  rigid 
sphere  and  deforming  cube  programs  and  the  backface  model  of 
Eqs .  (44)  and  (45).  In  all  cases,  6=4.  The  theoretical 
values  of  Ej.  ,  which  have  been  verified  for  each  material  in  the 

iv 

semi-infinite  case,  were  used  in  this  calculation.  In  the  same 
figure  are  plotted  the  predicted  residual  velocities  for  no  back- 
face  effect.  The  value  of  target  thickness  at  which  the  resi¬ 
dual  velocity  equals  zero  measures  the  stopping  thickness  re¬ 
quired  for  each  projectile  at  the  indicated  initial  velocity. 
Arrows  on  each  figure  indicate  the  experimental  thickness  which 
stopped  or  failed  to  stop  the  projectile.  We  conclude  that  the 
backface  model  with  6=4  works  quite  well  for  RHA,  titanium 
and  5083  aluminum. 

It  should  be  pointed  out  that  there  are  exceptions  to  this 
model.  Certain  composite  materials,  such  as  fiberglass  and 
Kevlar  woven  rovings  are  better  modeled  as  having  no  backface 
effects,  or  6  <<  1  .  Similarly,  brittle  materials  act  as  though 
6  >>  4  .  Fortunately,  however,  a  large  number  of  ductile 
materials  including  many  important  armor  materials  are  described 
by  6  =  4. 
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V.  OBLIQUE  PENETRATION 

When  a  rod  impacts  an  armor  plate  at  an  oblique  angle,  the 
forces  on  the  rod  will  not  be  axially  symmetric.  Thus  we  include 
in  the  rod  program  a  lateral  force  F^  acting  on  the  head  of 
the  rod,  as  well  as  the  axial  force  ,  which  was  described  in 

Sections  II  and  III.  A  bending  mode,  characterized  by  Ug  ,  the 
lateral  displacement  of  the  head  relative  to  the  axis  of  the 
shaft,  and  a  twisting  angle  f  relative  to  the  direction  of  rod 
motion  are  included.  In  addition,  the  trajectory  of  the  rod  no 
longer  will  follow  the  initial  direction  of  flight,  so  instead 
of  one  parameter  p  for  penetration  we  employ  p  as  the  total 
length  of  penetration  plus  X  ,  the  angle  of  penetration  relative 
to  the  original  velocity  direction.  All  these  quantities  are  de¬ 
fined  in  Figs.  10a  and  10b. 

First,  consider  the  dynamics  of  the  rod  itself.  The  total 
energy  of  the  rod  is 

T  -  Ws+V(l  +  <46) 

where  I  is  the  rod  moment  of  inertia,  U  is  the  center  of 
mass  velocity  of  the  rod  and  Ug  the  bending  velocity.  Since 
bending  and  other  nonaxial  effects  are  generally  small  correc¬ 
tions  to  the  total  penetration,  we  treat  the  rod  as  a  single 
element,  not  separating  it  into  head  and  shaft,  for  the  purpose 
of  calculating  these  effects. 

The  variables  expressing  lateral  deflection  are  then 
determined  by 

'i  =  Fl  (47) 


where  L  is  the  total  length  of  the  rod,  and 

i  -  5FL 
B  M  +  M- 
s  r 


(48) 


The  lateral  force  Fg  is  computed  in  the  following  way. 
As  the  rod  impacts  a  target  at  an  oblique  angle,  one  comer  of 


U0)  initial  velocity 
U,  instantaneous 
velocity 


0 ,  bending  dis¬ 
placement 


the  front  end  first  makes  contact  with  the  target  when  the 
center  of  the  front  face  of  the  rod  is  still  a  distance 


d  =  a  tan  6 

from  the  target,  measured  in  the  direction  of  the  rod  velocity. 

We  assume  the  plastic  and  drag  forces  acting  on  the  rod  begin 
to  increase  from  zero  at  this  initial  instant  of  contact  and 
rise  linearly  to  their  fully- imbedded  values  when  the  face  of 
the  rod  is  fully  imbedded  in  the  target.  The  penetration  p  is 
measured  in  the  program  from  this  initial  contact  point.  At 
this  point  the  center  of  the  front  face  of  the  rod  is  a  distance 

(a  •  sin  6)  from  the  target,  measured  along  a  normal  to  the  target 
face.  Thus,  we  treat  the  target  as  having  an  effective  thickness 
t  =  T  +  a  sin  6  ,  where  T  is  the  true  target  thickness.  The 
total  plastic  and  drag  pressure  exerted  on  the  head  of  the  rod 
is  computed  from  (1),  and  and  F^  are  found  by  multiplying 

this  pressure  by  the  front  face  area  and  lateral  area,  respec¬ 
tively,  in  contact  with  the  target.  The  head  of  the  rod  will 
begin  to  widen  as  soon  as  the  decelerating  force  of  the  target 
acts  on  it.  Full  embedding  of  the  head  (see  Fig.  10)  occurs 
when 

p  =  pQ  e  (a  +  b)  tan  (6  -  x) 

The  pressure  at  the  contact  interface  is 

Q 

Pt(E* t  +  ^<Vf  +  2  )2) 

and  according  to  our  assumption  of  a  linear  initial  increase  in 
the  contact  area,  the  force  exerted  on  the  front  face  of  the  rod 

is  C 

FA  “  ^  "b2pt(E*t  +  ^(Vf  +  O2),  p  <  po  (50) 

and  the  force  exerted  on  the  side  face  is 

fl  -  jp  4iTJlbpt(E*  +  ^(Vf  +  £) 2 )  ,  p<p0  (51) 
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We  assume  this  linear  increase  in  contact  area  holds  true  up  to 
the  point  of  full  embedding  at  the  penetration  p  =  pQ  .  For 
p  >  pQ  ,  F^  drops  to  zero  since  the  lateral  forces  acting  on 
all  sides  of  the  head  cancel  each  other,  yielding 


FA  =  *b2pt(E*t  +  ^(Vf  +  SI)2) 


FL  =  0 


(52) 


Upon  exiting  from  the  backface  of  the  target,  the  backface  model 
described  in  the  previous  section  is  generalized  to  oblique  exit 
in  the  following  way.  E,vt  and  the  drag  coefficient  near 

the  backface  are  assumed  to  decrease  according  to  (45) ,  where  the 
distance  from  the  penetrator  to  the  backface  is  taken  as  the  pro¬ 
jected  distance  measured  normal  to  the  backface.  The  reduced 
values  of  E.v  and  CD  are  substituted  directly  into  the 
formula  for  FA  .  Hear  the  backface,  the  projected  distance  from 
the  center  of  the  rod  face  to  the  backface  is  x  -  p  cos  (8  -  x)  » 
and  the  expression  for  the  effective  corresponding  to  (45) 

E*t  -  <*E*  (53) 


where  a  is  defined  by 

1  ,  x  -  p  cos  (8  -  x)  >  Sr 

a  '  r  ,  T  ■  p  C0s(e  -  X)  <  6r  (54) 

A  similar  dependence  is  assumed  for  near  the  backface. 

x  is  the  effective  plate  thickness  (T  +  a  sin  8)  measured  from 
the  point  of  initial  contact.  As  the  axial  force  decreases  near 
the  backface,  the  lateral  force  F^  increases  because  of  the 
imbalance  in  the  effective  E*  of  the  target  on  the  sides  of  the 
head.  Therefore,  we  assume 

FL  «  <1  -  a)(4*b£pt(E*t  +  ^2(Vf  +  £)2))  (55) 

These  assumptions  completely  specify  the  backface  effects  for 
oblique  exit  from  the  target. 


The  equations  above,  in  addition  to  those  described  in 
Section  III,  are  integrated  numerically  to  predict  the  residual 
mass,  residual  velocity,  ballistic  limit  velocity,  and  other 
parameters  for  oblique  penetration  as  well  as  normal  penetration 
The  oblique  model  described  here  has  been  included  in  the  PE. '5 
code,  but  has  not  thus  far  been  added  to  the  code  ROD.  A  copy 
of  the  code  PER  is  reproduced  in  Appendix  II.  The  input  para¬ 
meters  vrtiich  require  specification  are  the  geometric  dimensions, 
densities  and  E,v  values  of  the  target  and  penetrator  materials 
and  the  initial  striking  velocity  and  olbiquity  angle  of  the 
penetrator.  The  output  includes  ballistic  limit  velocity  or 
penetration  depth,  residual  mass  and  residual  velocity  of  the 
penetrator . 
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VI .  COMPARISON  OF  EXPERIMENT  WITH  THEORY 


In  order  to  check  the  programs  ROD  and  PEN  against  experi¬ 
ment  under  controlled  conditions,  several  long  rods  of 

L/D  =  10  were  fired  for  us  into  Rolled  Homogeneous  Armor 

2 

targets  of  Brinell  hardness  290  kg/mm  by  the  Ballistics 
Research  Lab,  Aberdeen  Proving  Grounds.  In  order  to  avoid  back- 
face  effects  which  might  modify  the  effective  E.,,  of  the  target 
in  these  initial  experiments,  very  thick  targets  of  thickness 
greater  than  twice  the  total  rod  penetration  were  used.  The  rods 
were  chosen  to  provide  a  variety  of  materials  and  strengths,  from 
1018  steel  to  soft  lead  to  Mallory  3000,  a  tungsten  alloy.  The 
resulting  data  were  compared  to  predictions  from  the  ROD  and  PEN 
programs  to  select  a  best  fit  value  of  c  ,  which  characterizes 
the  maximum  head  width,  and  x  .  which  relates  E,tj  to  the  rod 
strength  in  Eq .  131).  The  theoretical  values  of  E,.,  for  the  tar¬ 
get  and  penetrator  material  were  used,  based  on  formula  (30)  and 
the  melting  temperature,  hardness  and  heat  capacity  of  the 
respective  materials.  The  best  fit  values  for  ROD 
were  found  to  be  £  *  1.36  ,  >  =  .42  ,  and  for  PEN,  c  =1.7 

and  x  =  1 ■ 0  • 

A  comparison  of  the  resultant  theoretical  predictions  with 
experiment  for  the  code  ROD  is  shown  in  Figs.  11  through  13.  The 
corresponding  fit  to  the  data  for  the  code  PEN  is  shown  in  Fig.  14. 
The  high  velocity  lead  rod  deformed  upon  exit  from  the  gun  barrel 
and  had  a  highly  irregular  shape  and  L/D  ~  5  upon  impact  at  the 
target.  We  have  used  L/D  =5  in  computing  the  theoretical 
penetration  for  this  data  point.  Agreement  with  experiment  in  all 
cases  is  within  about  1570. 

Next,  the  code  was  tested  for  rods  against  finite  thickness 
targets  at  normal  incidence.  The  values  of  the  parameters  x 
and  £q  found  above  were  kept  the  same.  In  the  numerical  code 
this  set  of  experiments  amounts  to  a  test  of  the  accuracy  of  the 
backface  model,  presented  in  Section  IV.  The  ballistic  limit 
velocity  Vg^  was  determined  in  the  code  by  incrementally 
raising  the  striking  velocity  of  the  rod  until  penetration  was 
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achieved.  This  theoretical  Vt 


from  the  ROD  code  was  cotnoared 


to  one  set  of  experiments  (Lambert^)  with  65  gm  Bearcat  steel 

rods  of  L/D  =  5  ,  10  and  20  into  RHA  targets  of  various  hard- 

2 

ness  from  BHN  =  260  to  375  kg/mm  These  date  are  presented  in 
Figs.  15  through  17.  Agreement  is  excellent,  although  there  is 
a  slight  tendency  for  the  higher  L/D  rod  to  overpenetrate  and 
the  lower  L/D  rod  to  underpenetrate.  In  Figs.  18  through  20 
these  data  of  Lambert  are  plotted  in  dimensionless  form  against 
another  set  of  experiments  (Herr^)  in  which  1.94  ,  3.89  ,  and 
7.78  gram  Bearcat  steel  rods  were  fired  into  RHA  plates  which 
were  annealed  to  a  Brinell  hardness  of  400.  In  Figs.  18  through  20 
the  plate  thickness  is  expressed  in  units  of  T/D  or  (thickness)/ 
(rod  diameter) ,  since  it  is  the  rod  width  D  that  sets  the  dimen¬ 
sion  for  backface  effects.  It  is  apparent  in  the  figures  that  the 
two  sets  of  experimental  data  do  not  overlap,  particularly  for 
L/D  =  10  and  20  .  This  can  be  shown  to  be  a  result  of  the  dif¬ 


ferent  Brinell  hardnesses  of  the  armor  targets  used  in  the  two 
sets  of  experiments.^  Note,  for  example,  that  the  1.94  ,  3.89 
and  7.78  gram  rods  do  fall  on  the  same  ballistic  curve.  In 
Herr's  experiments,  all  targets  were  heat  treated  to  a  uniform 
hardness  of  BHN  =  400  .  Thus,  they  all  had  the  same  value  of 
Ej.  ,  about  215  Btu/lb.  In  Lambert's  data,  the  Brinell  hardness 
varies  from  about  260  to  375  km/mm^  with  corresponding  values 
of  E...  ranging  from  170  to  210  Btu/lb.  The  theoretical  curves 
corresponding  to  these  two  sets  of  E,v  values  are  plotted  in 
Figs.  18  through  20,  where  the  average  E*  of  Lambert's  targets, 
190  Btu/lb,  is  used.  The  theory  clearly  shows  the  same  shift 
in  VgL  with  hardness  that  is  seen  in  the  experimental  data. 
Thus,  the  different  values  of  E*  corresponding  to  different 
Brinell  hardnesses,  as  predicted  by  Eq .  (30),  account  for  the 
apparent  failure  of  scaling  in  the  data.  This  set  of  experiments 
provides  a  sensitive  test  of  the  ability  of  the  Integral  Theory 
to  predict  the  effect  of  materials  properties,  such  as  hardness, 
on  penetration  and  ballistic  limit. 
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In  Figs.  21  through  24  the  theoretical  value  of  the  residual 

9  10 

mass  versus  striking  velocity  is  compared  to  experiment  ’ 

for  Bearcat  steel  and  Mallory  3000  rods.  The  residual  mass 
computed  in  the  program  is  the  sum  of  the  mass  of  the  head  and 
shaft  of  the  rod  at  the  instant  the  target  backface  is  reached. 
Although  there  is  a  fair  amount  of  scatter  in  the  data,  good 
qualitative  agreement  is  attained  over  the  range  of  velocities 
and  target  thicknesses  in  the  experiment.  In  Figs.  25  through 
29  the  predicted  residual  velocity  of  the  rod  fragment  is  com¬ 
pared  with  experiment  for  various  thickness  targets  and  rods 
of  various  L/D  ,  for  Mallory  3000  as  well  as  steel.  These 
comparisons  are  typical.  The  agreement  is  good  between  the 
code  and  experiment  to  about  ±15%. 
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In  Figs.  30  and  31  typical  oblique  rod  shots  are  compared 
to  the  PEN  code  prediction.  The  ballistic  limit  is  plotted 
versus  striking  velocity  for  Bearcat  rods  into  RHA  at  60° 
incidence.  Again,  the  agreement  between  theory  and  experiment 
is  good  to  better  than  10%. 

In  order  to  provide  further  confirmation  of  the  numerical 
codes,  especially  over  a  wider  range  of  materials,  published 
data  were  obtained  for  long  rods  (wires)  of  gold,  tin,  aluminum, 
and  magnesium  fired  into  7075-T6  Aluminum  semi-infinite  targets. ^ 
The  values  of  E,v  for  the  targets  and  penetrators  were  obtained 
from  handbook  data  on  the  materials  involved,  and  substituted 
into  formula  (30)  .  The  resulting  values  of  E.,,  are  displayed 
in  Fig.  32,  together  with  the  experimental  and  theoretical 
curves  of  penetration  versus  velocity  for  the  four  rod  materials . 
The  penetration  is  normalized  to  allow  presentation  of  all 
curves  on  the  same  graph.  The  very  good  agreement  over  a 
range  of  rod  density  from  1.8  to  19  gm/cc  confirms  that  density 
variations  are  taken  into  account  correctly  in  the  code. 

As  yet  another  test  of  the  influence  of  the  material 
hardness  on  penetration,  published  data^  for  steel  and  Densimet 
17  rods  fired  into  various  types  of  semi-infinite  steel  targets 
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were  compared  to  the  PEN  numerical  code.  For  these  data,  the 
maximum  velocity  is  3000  m/sec,  considerably  higher  than  for 
the  previous  data  for  steel  rods  in  Fig.  11.  As  Figs.  33,34  show 
agreement  is  good  over  the  entire  range  of  velocity  for  all  of 
the  materials. 

Finally,  the  rod  code  was  tested  for  the  interesting  case 

of  a  dense,  strong  rod  into  a  low  density  target,  in  this  case 

12 

Bearcat  steel  rods  into  a  6061-T6  A1  target.  The  steel  rod, 
if  sufficiently  strong,  will  not  deform  at  low  velocities  and 
will  follow  the  ballistic  curve  of  a  nondeforming  projectile. 

At  high  enough  striking  velocity,  the  front  face  pressure  on  the 
rod  exceeds  the  yield  strength,  and  the  rod  begins  to  erode 
rapidly  and  behaves  like  a  deforming  rod.  In  Fig.  35,  data 
are  presented  which  display  this  transition  for  rods  of  various 
strengths.  The  theoretical  curves  corresponding  to  rods  of 
Rockwell  hardness  R  =  40  ,  50  and  60  (E*  =  207,  236  and  256 
Btu/lb)  clearly  show  that  the  stronger  rod  begins  to  deform 
at  higher  velocity.  These  data  provide  a  critical  test  of  the 
way  the  rod  strength  enters  the  code,  since  the  transition  from 
nondeforming  to  deforming  rod  is  so  sharp.  Experiments  such 
as  this  one  provide  a  sensitive  means  for  empirically  determining 
the  uniaxial  adiabatic  yield  strength  of  a  rod. 


VII.  CONCLUSIONS 


The  Integral  Theory  of  Impact  has  been  applied  to  the  pro¬ 
blem  of  modeling  the  behavior  of  long-rod  penetrators.  A  two¬ 
cell  model  for  the  deforming  rod  is  employed,  assuming  conser¬ 
vation  of  energy  and  momentum,  and  using  the  A.R.A.P.  concept  of 
adiabatic  hardness,  pE*  ,  to  account  for  material  strengths.  A 
numerical  code  based  on  these  assumptions  has  been  developed 
which  can  predict  the  performance  of  rod  penetrators  impacting 
finite  thickness  targets.  The  input  parameters  required  to 
operate  the  code  consist  only  of  the  physical  dimensions  of  the 
target  and  penetrator,  and  standard  handbook  properties  of  the 
materials,  such  as  density,  heat  capacity  and  Brinell  hardness, 
which  is  needed  to  compute  pE^  . 

The  code  predictions  of  ballistic  limit  velocity,  penetration 
depth,  residual  mass  and  residual  velocity  are  in  good  agreement 
(typically  ±15?0)  with  experiment  over  a  wide  range  of  materials 
and  velocities.  The  code  accurately  predicts  the  relative  im¬ 
provement  in  performance  of  a  rod  when  its  strength  is  increased, 
or  when  the  target  hardness  is  changed,  and  also  predicts  the 
approximate  velocity  at  which  a  rod  transitions  from  nondeforming 
deforming  penetration. 

The  current  treatment  of  oblique  impact  does  not  attempt  to 
handle  fracture  of  the  rod  shaft  or  jetting  of  the  rod  front  end 
during  impact.  However,  where  shaft  fracture  is  not  a  problem, 
the  code  predictions  of  ballistic  limit  are  in  good  agreement 
with  experiment . 
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APPENDIX  1 

LISTING  OF  ROD  CODE 
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